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Laminar  burning  velocities  and  flammability  limits  of  premixed  methane/air  flames  in  the  presence  of 
various  diluents  were  investigated  by  combined  use  of  experiments  and  numerical  simulations.  The 
experiments  used  a  1-m  free-fall  spherical  combustion  chamber  to  eliminate  the  effect  of  buoyancy, 
enabling  accurate  measurements  of  near-limit  burning  velocities  and  flammability  limits.  Burning  veloc¬ 
ities  were  measured  for  CH4/air  flames  with  varying  concentrations  of  He,  Ar,  N2  and  CO2  at  NTP.  The  lim¬ 
iting  concentration  of  each  diluent  was  measured  by  systematically  varying  the  composition  and  ignition 
energy  and  finding  the  limiting  condition  through  successive  experiment  trials.  The  corresponding  freely- 
propagating,  planar  1-D  flames  were  simulated  using  PREMIX.  The  transient  spherically-expanding 
flames  were  simulated  using  the  1-D  Spherical  Flame  &  Reactor  Module  of  COSILAB  considering  detailed 
radiation  models.  The  results  show  that  helium  exhibits  more  complex  limit  behavior  than  the  other  dil¬ 
uents  due  to  the  large  Lewis  number  of  helium  mixtures.  The  near-limit  helium-diluted  flames  require 
much  higher  ignition  energy  than  the  other  flames.  In  addition,  for  the  spherically  expanding  helium- 
diluted  flames  studied  here  (Le  >  1 ),  stretch  suppresses  flame  propagation  and  may  cause  flame  extinc¬ 
tion.  For  the  C02-diluted  flames,  the  flame  speed  predicted  by  the  optically-thick  model  based  on  the  Dis¬ 
crete  Transfer  Method  (DTW)  and  a  modified  wide  band  model  has  better  agreement  with  measurements 
in  the  near-limit  region.  A  significant  amount  of  heat  is  absorbed  by  the  dilution  gas  CO2,  resulting  in  ele¬ 
vation  of  temperature  of  the  ambient  gases.  The  optically-thick  model,  however,  still  overpredicts  flame 
speed,  indicating  a  more  sophisticated  radiation  property  model  may  be  needed.  Finally,  the  chemical 
effect  of  CO2  on  flame  suppression  was  quantified  by  a  numerical  analysis.  The  results  show  that  the 
chemical  effect  of  CO2  is  more  important  than  the  other  diluents  due  to  its  active  participation  in  the 
reaction  CO2  +  H  =  CO  +  OH,  which  competes  for  H  radicals  with  the  chain-branching  reactions  and  thus 
reduces  flame  speed. 

©  2010  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

The  inert  gases  such  as  helium,  argon,  nitrogen,  and  carbon 
dioxide  are  an  important  class  of  fire  extinguishing  agents  because 
they  are  natural,  non-ozone-depleting,  non-toxic  and  non-pyro- 
lytic  gases.  Such  gases  may  be  essential,  for  example,  for  extin¬ 
guishing  fires  in  confined  environments  such  as  spacecraft.  The 
extinguishing  mechanisms  of  these  diluents  have  been  studied 
experimentally  and  computationally  for  various  flame  configura¬ 
tions  at  terrestrial  or  microgravity  conditions  [1-8]. 

Laminar  burning  velocity  and  flammability  limits  of  premixed 
flames  are  key  parameters  of  both  fundamental  and  practical  sig¬ 


*  Corresponding  author.  Address:  School  of  Aeronautics  &  Astronautics,  Purdue 
University,  701W.  Stadium  Ave.,  West  Lafayette,  IN  47907-2045,  USA.  Fax:  +1  765 
494  0307. 

E-mail  address:  lqiao@ecn.purdue.edu  (L.  Qiao). 


nificance.  Accurate  measurement  of  unstretched  laminar  burning 
velocity  is  necessary  for  assessing  combustion  theories  and  the  val¬ 
idation  of  numerical  models.  Numerous  studies  have  shown  that 
gravity  has  significant  influence  on  near-limit  behavior  of  laminar 
and  turbulent  flames  [9-13].  Ronney  [13]  provided  a  comprehen¬ 
sive  review  on  premixed-gas  flames  in  microgravity.  Generally, 
the  flammability  limits  obtained  at  1-g  and  p-g  are  different.  Burn¬ 
ing  velocities  of  near-limit  flames  are  difficult  to  measure  at  1-g 
because  the  natural  convection  causes  severe  flame-front  distor¬ 
tion.  For  example,  the  outwardly-propagating  spherical  flames, 
which  have  been  traditionally  used  to  measure  flame  speed,  would 
propagate  upward  at  limit  conditions  and  form  a  mushroom  shape 
due  to  buoyant  convection. 

In  a  previous  study  [14]  we  presented  results  from  a  1-m  free- 
fall  experiment  that  eliminated  the  effect  of  buoyancy  on  premixed 
flames.  Those  experiments  used  a  stationary  shadowgraph  system 
to  image  flame  propagation  inside  a  dropping  spherical  chamber. 
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Nomenclature 

D 

mass  diffusivity 

a 

Stefan-Boltzmann  constant 

K 

flame  stretch 

kp 

Planck  mean  absorption  coefficient 

L 

Markstein  length 

V 

spectrum  wave  number 

Ma 

Markstein  number,  LISq 

M 

total  number  of  bands 

characteristic  flame  thickness,  DuISloo 

N 

total  number  of  rays 

rf 

flame  radius 

n 

grid  cell 

t 

time 

kv 

spectral  absorption  coefficient 

Sl 

local  laminar  burning  velocity 

KV 

optical  depth 

Sloo 

planar  unstretched  laminar  burning  velocity 

hv 

spectral  blackbody  radiation  intensity 

P 

pressure 

rj 

the  spectrum 

T 

temperature 

Qr 

radiation  heat  source 

To 

ambient  temperature 

I 

spectral  blackbody  radiation  intensity 

X 

mole  fraction 

Q 

solid  angle 

P 

density 

^n,ij 

radiation  intensity  of  band  j  and  ray  i  at  grid  n 

fuel-equivalence  ratio 

The  near-limit  flames  maintained  a  spherical  shape  that  allowed 
their  laminar  burning  velocity  to  be  inferred  from  the  measured 
flame  radius  versus  time.  Measurements  of  near-limit  burning 
velocities  under  microgravity  conditions  for  H2/air/diluent  pre¬ 
mixed  flames  at  various  fuel-equivalence  ratios  and  pressures 

[15]  were  reported  using  this  facility.  The  measured  velocities 
were  compared  to  calculations,  which,  when  analyzed,  revealed 
that  inaccuracies  in  the  three-body  termination  rates  for 
H  +  O2  +  M  =  HO2  +  M  reactions  and  in  mass  diffusion  coefficients 
for  H2  diffusion  are  the  most  likely  explanation  for  the  near-limit 
differences.  In  these  studies  [14,15],  numerical  simulations  were 
performed  for  1-D  planar  freely-propagating  flames  using  PREMIX 

[16] .  Radiative  heat  loss  was  either  not  considered  or  was  simu¬ 
lated  using  a  simple  optically-thin  model.  The  results  show  that 
radiative  heat  loss  was  insigniflcant  for  He,  Ar,  or  N2-diluted  hydro¬ 
gen/air  flames.  It,  however,  could  be  signiflcant  for  C02-diluted 
flames.  The  reabsorption  effect  of  CO2  can  affect  near-limit  behav¬ 
ior  when  CO2  concentration  is  high.  Furthermore,  the  1-D  planar 
flame  simulation  using  PREMIX  may  not  represent  the  spherical 
flames  as  investigated  in  the  experiments,  which  are  subject  to 
flame  stretch.  Stretch,  combined  with  radiation/reabsorption, 
could  change  flame  dynamics  and  extinguishing  processes.  There¬ 
fore,  a  more  sophisticated  modeling  that  combines  all  these  effects 
is  desirable  for  a  better  interpretation  of  flame  extinction  behavior 
by  diluents. 

The  combined  effects  of  stretch  and  radiation  on  premixed 
flames  have  been  studied  by  several  researchers  [17-35].  Farmer 
and  Ronney  [17]  numerically  examined  the  dynamics  properties 
of  near-limit  CH4/air  flames  in  planar  and  spherical  geometries 
considering  radiative  heat  loss  from  CO2  and  H2O.  Their  results 
show  that  spherical  geometry  simulations  had  good  qualitative 
and  quantitative  agreement  with  results  obtained  in  microgravity 
experiments.  The  results  also  show  that  the  extinguishment  pro¬ 
cesses  at  microgravity  are  a  result  of  the  interactions  among  the 
Lewis  number  effect,  flame  front  curvature  and  heat  loss  due  to 
gas  radiation.  Buckmaster  [22]  examined  the  behavior  of  twin  pre¬ 
mixed  opposed-flow  flames  considering  radiation  heat  loss.  The 
numerical  simulations  captured  a  number  of  different  flame  re¬ 
sponses  for  a  wider  range  of  mixture  strength.  Bechtold  et  al. 
[32]  developed  a  general  theory  of  non-adiabatic  premixed  flames 
of  arbitrary  shape  that  fully  accounts  for  the  hydrodynamic  and 
diffusive-thermal  processes  and  incorporates  the  effects  of  volu¬ 
metric  heat  losses.  Wang  and  Niioka  [31]  numerically  investigated 
radiation  reabsorption  effects  in  CH4-air  counterflow  premixed 
flames  and  found  that  the  reabsorption  of  emitting  radiation  leads 
to  substantially  wider  flame  thickness  and  higher  flame  temperature 


than  those  calculated  by  using  the  optically-thin  model.  Ju  and  his 
co-workers  [23-30,34,35]  conducted  extensive  studies  on  the 
combined  effects  of  radiation  and  stretch  on  flammability  limits 
for  various  Lewis  number  mixtures.  In  particular,  their  recent  work 
[34]  show  that  the  flame  geometry  has  a  signiflcant  impact  on 
radiation  absorption  and  that  the  1-D  planar  radiation  model  was 
not  valid  for  the  computation  of  the  flame  speed  of  a  spherical  flame. 

To  validate  radiation  models  and  to  understand  the  dynamics 
and  extinction  process  of  premixed  flames,  experimental  data  such 
as  flame  speed  are  useful,  which,  however,  are  not  easy  to  obtain 
for  near-limit  flames  because  of  buoyancy.  In  the  present  paper, 
measurements  of  laminar  burning  velocities  and  extinction  limits 
for  methane/air/diluent  flames  using  the  same  1-m  free-fall  facility 
were  reported.  The  transient  outwardly-propagating  spherical 
flames  were  simulated  using  the  Spherical  Flame  Module  of  COSl- 
LAB  [36]  with  detailed  radiation  models.  The  objective  was  to 
understand  how  and  when  the  premixed  methane-air  flames  are 
extinguished  by  various  inert  gases  (He,  Ar,  N2  and  CO2),  and  in 
particular,  to  understand  the  effects  of  Lewis  number,  flame 
stretch,  and  radiation  emission  and  absorption  of  CO2  on  flame 
propagation  and  extinction.  Here  measurements  of  burning  veloc¬ 
ities  of  stoichiometric  methane/air/diluent  mixtures  at  NTP  at 
varying  diluent  concentrations  are  first  presented.  The  measured 
velocities  are  then  compared  with  numerical  simulations.  After 
that  the  measured  concentration  of  each  inert  gas  at  the  point  at 
which  the  flame  is  extinguished  is  reported  and  is  compared  with 
the  data  in  the  literature.  This  limiting  concentration  was  deter¬ 
mined  by  systematically  varying  the  composition  and  ignition  en¬ 
ergy  and  finding  the  limiting  condition  through  successive 
experiment  trials.  The  results  show  that  helium  behaves  in  a  com¬ 
plex  way  at  the  limit  due  to  its  high  diffusivities.  Then  the  effects  of 
radiation  emission  and  absorption  on  the  C02-diluted  flames  are 
discussed.  The  results  show  that  the  flame  speed  predicted  based 
on  the  Discrete  Transfer  Method  (DTM)  and  a  modified  wideband 
model  is  higher  than  that  predicted  by  the  optically-thin  limit 
model.  A  signiflcant  amount  of  radiation  heat  loss  from  the  hot 
products  is  reabsorbed  by  the  cold  ambient  gases,  resulting  in  tem¬ 
perature  increase  in  the  unburned  gas  and  the  pre-heat  zone, 
which  consequently  increases  reaction  rate  and  flame  speed.  Final¬ 
ly,  the  pure  chemical  effect  of  CO2  in  flame  suppression  was  quan¬ 
tified  by  a  numerical  analysis.  The  results  show  that  the  chemical 
effect  of  CO2  can  cause  as  much  as  30%  reduction  in  flame  speed 
near  extinction  limits.  This  is  because  the  competition  of  CO2  for 
H  radical  through  the  reaction  CO  +  OH  =  CO2  +  H  with  the  chain¬ 
branching  reaction  H  +  O2  =  0  +  OH  reduces  H,  OH,  and  0  radical 
concentrations,  leading  to  lower  reaction  rate  and  flame  speed. 
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2.  Experimental  method 


The  experimental  methods  are  similar  to  earlier  work  [14,15] 
and  will  be  described  only  briefly  here.  The  free-fall  facility  con¬ 
sists  of  a  support  tower,  a  free-falling  spherical  combustion  cham¬ 
ber,  a  spark  generator,  a  deceleration  box,  and  shadowgraph  optics 
that  record  the  flame  propagation  as  the  chamber  falls.  The  idea 
was  to  use  a  stationary  optical  system  to  image  flame  propagation 
inside  a  dropping  chamber.  The  free-fall  chamber  was  held  at  the 
top  of  the  tower  by  an  electromagnet  before  being  dropped.  As 
the  chamber  was  released,  a  hall-effect  sensor  detected  the  motion 
and  sent  a  trigger  pulse  to  the  delay  generator.  A  short  delay  was 
then  provided  to  allow  oscillations  from  the  chamber  release  to  de¬ 
cay  to  sufficiently  low  levels.  After  the  delay,  the  generator  sent  out 
two  signals,  one  of  which  triggered  a  high-voltage  spark  generator 
connected  to  electrodes  that  ignited  the  mixture,  and  the  other 
triggered  a  high-speed  digital  video  camera  which  recorded  shad¬ 
owgraph  images  of  the  flame  propagation  within  the  free-falling 
chamber.  The  spark  gap  and  spark  energy  were  adjusted  to  be  as 
close  as  possible  to  the  minimum  ignition  energies. 

The  present  measurements  were  restricted  to  flames  having 
diameters  10  mm  <  d  <  60  mm.  The  lower  limit  provides  sufficient 
time  for  disturbances  introduced  by  ignition  to  decay,  while  the 
upper  limit  limits  the  pressure  rise  in  the  chamber  during  the  mea¬ 
suring  period  to  less  than  0.7%  of  the  initial  pressure.  No  results 
were  considered  where  the  flame  surface  appeared  wrinkled  due 
to  flame  instabilities.  As  in  previous  experiments,  the  measure¬ 
ments  presented  here  were  limited  to  0.05  where  So  is  the 
characteristic  flame  thickness  and  r/is  flame  radius.  This  minimizes 
effects  of  curvature  and  transient  phenomena  associated  with  large 
flame  thickness  during  the  early  stages  of  flame  formation.  Under 
these  assumptions,  quasi-steady  expressions  for  the  local  laminar 
burning  velocity  and  flame  stretch  are  given  by  [37] 
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2  ^ 

Pu  dt  ’  rfdt’ 


(1) 


where  Si  is  the  observed  stretched  flame  propagation  speed  into  the 
unburned  gas  and  K  is  the  flame  stretch.  The  density  ratio  needed  to 
determine  Si  was  computed  assuming  adiabatic  constant-pressure 
combustion  with  the  same  concentrations  of  elements  in  the  un¬ 
burned  gases.  These  calculations  were  carried  out  using  the  adiabatic 
equilibrium  algorithms  of  McBride  et  al.  [38].  The  flame  propagation 
velocity,  drf/dt,  was  obtained  from  the  flame  radius  r/t)  as  measured 
from  the  shadowgraph  images  along  the  direction  perpendicular  to 
the  spark  electrodes,  where  disturbances  of  the  flame  surface  by 
the  electrodes  were  minimal.  From  Markstein  [39]  and  Clavin  [40], 
the  laminar  burning  velocity  Si  is  related  to  the  flame  stretch  K  for 
small  to  moderate  values  of  curvature  and  stretch. 


Si  =  Si^-LK.  (2) 

The  Markstein  length  L  is  a  measure  of  the  flame  response  to  the 
stretch  rate  K,  and  it  can  be  either  positive  (stable  to  preferential- 
diffusion)  or  negative  (unstable  to  preferential-diffusion),  depend¬ 
ing  on  the  reactants.  The  unstretched  laminar  burning  velocities 
Sloo  were  obtained  by  extrapolating  the  measured  stretched  lami¬ 
nar  burning  velocities  Si{I<)  as  a  function  of  the  local  stretch  rate 
K  to  zero  stretch.  The  Markstein  number  Ma  was  obtained  by  nor¬ 
malizing  L  with  a  characteristic  flame  thickness  So  deflned  as  DJ 
5loo.  where  Du  is  the  mass  diffusivity  of  fuel-oxidizer.  For  simplicity 
and  consistency,  we  used  the  binary  diffusivity  of  CH4-O2 
(21.9  mm^/s)  for  all  mixtures  including  the  helium-diluted  flames, 
even  though  helium  addition  can  change  the  mass  diffusivity  of  the 
mixture  significantly. 

Experimental  conditions  for  CH4/air/diluent  flames  are  sum¬ 
marized  in  Table  1,  which  shows  reactant  mixtures  at  room 


temperature  (293  ±  0.5  K),  a  fuel-equivalence  ratio  of  1.0,  a  pres¬ 
sure  of  1.0  atm,  and  diluent  concentrations  from  0%  to  the  extinc¬ 
tion  limit  for  He,  Ar,  N2  and  CO2. 


3.  Computational  method 

Numerical  simulations  of  the  corresponding  planar  unstretched 
premixed  flames  were  carried  out  using  the  steady,  1-D  laminar 
premixed  flame  code  PREMIX  [16].  CHEMKIN  is  used  as  a  prepro¬ 
cessor  to  determine  thermochemical  and  transport  properties  from 
the  database  of  Kee  et  al.  [41,42].  The  simulations  used  the  detailed 
chemical  reaction  mechanism  GRl-Mech  3.0  [43],  which  involves 
325  reactions  and  53  species.  Multi-component  diffusion  and  var¬ 
iable  thermodynamic  properties  were  considered.  The  resulting 
computed  unstretched  burning  velocities  were  compared  for  each 
diluent  type  and  concentration  with  the  corresponding  measured 
values  obtained  as  described  above. 

To  evaluate  the  combined  effects  of  chemical  kinetics,  radiation, 
and  stretch  on  methane/air/diluent  flames,  we  also  simulated  the 
outwardly-propagating  spherical  flames  using  the  1-D  Spherical 
Elames  Module  of  COSILAB,  whose  original  version  was  called 
RUN  IDE  and  was  developed  by  Rogg  and  coworkers  at  Cambridge 
University  [44].  The  simulations  are  for  reactant  mixtures  initially 
at  room  temperature  293  K  and  are  for  flame  propagation  at  con¬ 
stant  pressures.  Similar  to  the  1-D  planar  flame  simulations,  GRl- 
Mech  3.0  was  used  in  the  transient  simulations.  CHEMKIN  was 
used  as  a  preprocessor  to  determine  thermochemical  and  transport 
properties  from  the  database  of  Kee  et  al.  [16,41,42,45].  Eor  the 
effective  solution  of  time-dependent  reactive-flow  problems  an 
extrapolation  method  called  EULER  [36]  was  used,  which  is  based 
on  algorithms  suitable  for  the  solution  of  differential  algebraic  sys¬ 
tems.  An  adaptive  gridding  method  was  also  used  to  improve  com¬ 
putational  efficiency. 

The  computational  domain  was  set  to  have  a  radius  of  8  cm.  The 
reactant  mixture  was  numerically  ignited  by  a  hot  spot  located  at 
the  origin.  The  size  of  the  hot  spot  is  typically  below  1/10  of  the 
computational  domain,  with  larger  spots  used  for  near-extinction 
flames.  After  the  reactant  mixture  was  ignited,  a  propagating  flame 


Table  1 

CH4/air/diiuent  laminar  premixed  flame  test  conditions.^ 


Diluent 

Xd 

PulPb 

Sloo  (cm/s) 

Le^ 

- 

0.00 

7.51 

37.0 

0.99 

He 

0.10 

7.26 

32.0 

1.70 

He 

0.20 

6.92 

25.9 

2.40 

He 

0.30 

6.55 

20.7 

3.11 

He 

0.35 

6.32 

17.1 

3.46 

Ar 

0.10 

7.26 

28.3 

0.99 

Ar 

0.20 

6.92 

20.3 

0.98 

Ar 

0.30 

6.55 

14.3 

0.98 

Ar 

0.40 

6.07 

9.0 

0.97 

Ar 

0.45 

5.79 

6.1 

0.97 

Ar 

0.48 

5.60 

4.3 

0.97 

N2 

0.10 

7.11 

25.5 

0.99 

N2 

0.15 

6.86 

20.7 

0.98 

N2 

0.20 

6.60 

15.9 

0.98 

N2 

0.25 

6.33 

11.8 

0.98 

N2 

0.30 

6.07 

8.6 

0.98 

N2 

0.35 

5.78 

5.4 

0.98 

CO2 

0.05 

7.15 

19.2 

0.97 

CO2 

0.10 

6.82 

11.4 

0.94 

CO2 

0.15 

6.46 

9.0 

0.92 

CO2 

0.20 

6.10 

5.5 

0.89 

CO2 

0.22 

5.96 

4.4 

0.89 

^  At  P  =  1 .0  atm,  (/)  =  1 .0,  and  298  ±  0.5  K. 

^  Lewis  number  is  evaluated  based  on  the  thermal  diffusivity  of  the  reactant 
mixture  and  the  binary  mass  diffusivity  of  CH4-O2  (21.9  mm^/s). 
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was  established,  ahead  of  which  was  the  quiescent  unburned  mix¬ 
ture  and  behind  which  were  the  combustion  products.  The  hot  spot 
temperature  profile  can  be  adjusted  so  that  the  amount  of  energy 
supplied  to  the  reactant  mixture  can  be  varied  to  be  close  to  the 
minimum  ignition  energy,  similar  to  the  experiment.  For  boundary 
conditions,  the  mole  fraction  of  each  gas  has  zero-gradient  at  the 
origin  and  has  fixed  values  at  the  right  boundary.  To  determine 
flame  propagation  speed,  we  chose  the  location  where  the  temper¬ 
ature  is  400  K  to  be  flame  front.  The  propagation  speed  of  expand¬ 
ing  flame  surfaces,  however,  does  not  depend  on  the  specific  choice 
of  flame  position  for  thin  flames  according  to  a  previous  study  [36]. 
The  methodology  to  obtain  unstretched  laminar  burning  velocity 
and  Markstein  length  from  the  outwardly-propagating  spherical 
flame  simulations  were  similar  to  the  experiments,  following  Eqs. 
(1)  and  (2).  To  remove  local  disturbances  introduced  by  the  numer¬ 
ical  scheme,  the  flame  speed  data  were  smoothed  using  LOWESS 
scheme  [46]  with  a  window  size  of  0.5.  To  avoid  ignition  distur¬ 
bances,  flame  speed  data  were  processed  for  r>  1.5  cm. 

Radiation  becomes  increasingly  important  when  flammability 
limits  are  approached.  Therefore  it  needs  special  consideration  in 
the  present  study  which  concerns  flame  extinction.  A  radiation 
model  typically  consists  of  a  radiation  property  model  for  the  radi¬ 
ative  properties  of  gases  and  a  radiative  transfer  model  which 
solves  the  radiative  transfer  equation  (RTE).  Here,  we  employed 
two  radiation  models  in  the  transient  simulations  using  COSILAB: 
an  opticahy-thin  limit  model  and  an  opticahy-thick  model  based 
on  the  Discrete  Transfer  Method  (DTM)  [47],  both  in  1-D  spherical 
coordinates.  The  optically-thin  limit  model  basically  assumes  that 
radiation  can  pass  through  the  medium  without  significant  absorp¬ 
tion.  With  gray  gas  assumption,  for  which  a  mean  absorption  coef¬ 
ficient  is  used  over  the  entire  wavelength  spectrum,  the  opticahy- 
thin  limit  model  can  be  simplified  as  [47] 

V-qr  =  k,{T)i4aT^-4aTl)  (3) 

where  V  Qr  is  the  radiative  heat  loss  rate  per  unit  volume,  kp  is  the 
Planck  mean  absorption  coefficient,  a  is  the  Stefan-Boltzmann  con¬ 
stant.  Radiation  of  the  gaseous  species  CO2,  H2O,  and  CO  were  con¬ 
sidered  and  their  absorption  coefficients  were  taken  from  [48]. 

The  details  of  the  Discrete  Transfer  Method  (DTM)  are  discussed 
in  Ref.  [47]  and  will  be  briefly  described  here.  The  DTM  is  based  on 
the  solution  of  selected  beams  of  radiation  and  on  the  subsequent 
computation  of  the  radiation  sources,  which  combine  the  finite  dif¬ 
ference  control  volumes  of  the  flow  procedure  due  to  the  passage 
of  the  beams.  The  individual  directions  of  rays  are  specified  in  ad¬ 
vance  rather  than  being  chosen  at  random.  The  change  in  spectral 
intensity  of  the  ray  passing  through  the  absorbing  and  emitting 
medium,  /v,  can  be  written  as 


dKy 


-Iv-hl 


bv 


(4) 


where  Ky  is  the  optical  depth  defined  as  Ky  =  kyds*  where  s*  =  0  is 
at  a  cell  entry  and  ky  denotes  the  spectral  absorption  coefficient  of 
the  gas  mixture,  v  is  spectrum  wavenumber  and  hy  is  the  spectral 
blackbody  intensity.  The  net  gain  of  radiation  energy  qr,n-ti2 
through  the  computational  cell  [r^.i,  r^],  whose  gradient  constitu¬ 
ent  a  source  term  in  the  energy  conservation  equation,  can  be  eval¬ 
uated  once  the  intensities  along  all  the  left  and  right  traveling  rays 
are  known.  Eor  spherical  geometry,  the  gradient  of  qr,n-i/2  is  ob¬ 
tained  by  integration  of  the  net  intensity  over  ah  solid  angles  in  a 
discretised  form  and  can  be  written  as  [47] 


(V  ■  Qr)n-t/2 


n  d(r^Qr}\ 

Vr2  dr  ) 

iCijrl  -  )  cos  ejAvjAQi 

(r„-r„_i)(r„  +  r„_i)^ 


(5) 


where  n  =  2, . . . ,  NG  (NG  is  the  number  of  grids),  In,ij  is  the  radiation 
intensity  of  ray  i  at  grid  n.  A6  is  given  by  (7r/2)/N,  0/  by  6i  =  (7r/2)/N-i 
and  AQi  =  sinOiAOi.  M  is  the  total  number  of  bands  and  N  is  the 
number  of  rays.  Eor  the  radiation  property  model  of  gas  species,  a 
modified  wide  band  model  was  used,  which  assumes  the  absorption 
and  emission  of  infrared  radiation  by  a  particular  gas  mainly  occur 
in  several  relatively  wideband.  Thirteen  bands  were  considered  for 
CO2,  CO  and  H2O  altogether.  They  are  the  2.0,  2.7,  4.3,  9.4, 10.4  and 
15  pm  vibrational-rotational  bands  of  CO2,  the  4.7  and  2.35  pm 
vibrational-rotational  bands  of  CO,  and  the  1.38,  1.87,  2.7  and 
6.3  pm  vibrational-rotational  bands  and  a  pure  rotational  band  of 
H2O. 


4.  Results  and  discussion 

4.1.  Laminar  burning  velocities 

Experiments  in  microgravity  are  essential  for  measuring  near- 
limit  burning  velocities.  Eor  slow-burning  mixtures  [14,15],  flames 
in  1-g  become  observably  buoyant,  with  the  flame  forming  a 
mushroom  shape  and  reaching  the  top  of  the  chamber  before  the 
bottom.  In  microgravity,  these  flames  remained  spherical  during 
propagation,  making  it  possible  to  measure  flame  radius  and  thus 
to  determine  flame  speed  and  stretch  rate.  Eor  ah  flames  consid¬ 
ered  here,  the  flame  surface  remained  smooth  at  the  early  stage 
of  the  propagation  process  and  no  preferential-diffusion  instability 
was  observed.  Ignition  sparks  might  result  in  irregular  ignition  ker¬ 
nels,  but  irregular  structures  quickly  smoothed  out  and  the  flame 
kernel  quickly  became  spherical.  Eor  very  slow  flames,  however, 
the  effect  of  irregular  ignition  kernels  was  larger  because  it  took 
a  longer  time  for  the  slow-propagating  surface  to  catch  up  with 
the  rest  of  the  surface. 

Eig.  1  shows  the  measured  and  computed  burning  velocities  as 
functions  of  the  molar  concentration  of  diluents  for  stoichiometric 
CH4/air  flames  at  normal  temperature  and  pressure  (NTP).  The  dil¬ 
uents  considered  are  N2,  Ar,  He,  and  CO2.  The  laminar  burning 
velocities  presented  here  are  unstretched  values,  5loo»  determined 
by  applying  Eq.  (2)  to  the  measured  values  of  Si  and  /C,  and  are  re¬ 
ported  in  Table  1.  The  estimated  Lewis  number  (Le)  is  also  reported 
in  Table  1  for  each  reactant  mixture.  The  computed  unstretched 
burning  velocities  were  obtained  from  both  the  steady-state  planar 
flame  simulations  using  PREMIX  [16]  and  the  spherically-expand- 


Fig.  1.  Measured  (symbols)  and  computed  (lines)  laminar  burning  velocities  as 
functions  of  the  mole  fractions  of  helium,  argon,  nitrogen  and  carbon  dioxide  for 
adiabatic,  premixed  CH^air/diluent  flames  at  ^  =  1.0  and  NTP.  Solid  lines:  simu¬ 
lation  results  for  the  planar  premixed  flames  using  PREMIX;  dashed  lines: 
simulation  results  for  the  spherically-expanding  flames  using  COSILAB. 
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ing  flame  simulations  using  COSILAB  [36].  The  computed  velocities 
shown  in  Fig.  1  are  for  adiabatic  flames  without  consideration  of 
radiation.  Effects  of  radiation  emission  and  absorption  on  flame 
properties,  particularly  for  near-limit  flames,  will  be  discussed  la¬ 
ter  in  Section  3.  For  the  COSILAB  simulations,  the  methodology  to 
obtain  Sloo  was  similar  to  the  experiment:  Si^o  was  obtained  by 
extrapolating  the  stretched  velocities  Sl  as  a  function  of  the  local 
stretch  rate  to  zero  stretch  (Eq.  (1)).  In  addition,  the  data  for  pro¬ 
cessing  were  limited  to  r>  1.5  cm  to  avoid  ignition  disturbances. 
The  density  ratios  needed  to  determine  Si  were  the  same  as  those 
for  the  experiments,  which  were  computed  assuming  adiabatic 
constant-pressure  combustion  using  the  adiabatic  equilibrium 
algorithms  of  McBride  et  al.  [38]. 

The  results  in  Fig.  1  show  that  all  diluents  considered  here  cause 
flame  speed  to  decrease  as  their  concentration  is  increased.  Similar 
to  the  conclusions  from  previous  studies  [3,14,15],  the  capabilities 
of  the  diluents  to  reduce  burning  velocity  at  any  flxed  diluent  con¬ 
centration  increases  in  the  order  He  <  Ar  <  N2  <  CO2,  with  CO2 
being  most  effective.  The  computed  Sloo  from  PREMIX  and  COSILAB 
simulations  generally  agree  well  with  each  other.  The  flame  speeds 
predicted  using  COSILAB  are  slightly  higher  than  those  using  PRE¬ 
MIX  for  the  helium  and  argon-diluted  flames. 

The  lowest  measured  burning  velocity  5loo  shown  in  Fig  1  is 
around  3-5  cm/s  for  the  CO2,  N2  or  Ar  diluted  flames  (except  He). 
Note  the  mole  fraction  of  the  diluent  that  corresponds  to  the  low¬ 
est  burning  velocity,  however,  is  not  the  measured  flammability 
limit,  which  is  a  topic  that  will  be  discussed  in  the  next  section.  He¬ 
lium  is  an  exception:  the  lowest  burning  velocity  we  could  mea¬ 
sure  for  the  helium-diluted  flames  is  about  15  cm/s  at  a  helium 
mole  fraction  of  38%,  above  which  the  ignition  kernel  dissipated 
and  the  flame  did  not  propagate  to  the  entire  chamber,  even 
though  much  higher  ignition  energy  was  used. 

4.2.  Markstein  lengths  and  numbers 

The  measured  and  computed  Markstein  length  (L)  and  Mark¬ 
stein  number  (Ma)  are  plotted  as  functions  of  the  mole  fractions 
of  diluents  for  stoichiometric  CH4/air/diluent  flames  in  Figs.  2 
and  3,  respectively.  As  shown  in  Table  1,  the  estimated  Lewis  num¬ 
ber  remains  almost  constant  with  N2-dilution  or  Ar-dilution,  and 
decreases  only  slightly  with  C02-dilution.  It,  however,  increases 
with  increasing  mole  fraction  of  helium. 

The  measured  Markstein  length  and  Markstein  number  have 
reasonable  agreement  with  the  computed  values  using  COSILAB. 
For  the  Ar,  N2,  and  C02-diluted  flames,  the  Markstein  length  varies 
little  with  increasing  mole  fraction  of  the  diluents,  except  that  a 
sudden  increase  was  observed  at  30-40%  N2-diltuion  near  extinc¬ 
tion  limit  (Fig.  2).  The  nearly  constant  Markstein  length  can  be  ex¬ 
plained  by  the  fact  that  addition  of  Ar,  N2,  and  CO2  does  not 
signiflcantly  alter  the  preferential  diffusion  behavior  of  the  mix¬ 
ture.  The  corresponding  Markstein  number  of  these  mixtures,  as 
shown  in  Fig.  3,  decreases  when  the  diluent  mole  fraction  in¬ 
creases.  This  is  because  the  flames  become  thicker  and  slower  with 
higher  diluent  concentrations.  The  trend  for  the  helium-diluted 
flames,  however,  is  different.  Both  the  Markstein  length  and  Mark¬ 
stein  number  increase  with  increasing  mole  fraction  of  helium.  The 
different  trend  observed  in  Figs.  2  and  3  is  attributed  to  the  much 
larger  diffusivity  of  helium  compared  to  O2,  N2,  Ar,  or  CO2. 

4.3.  Extinction  limits 

Research  on  flammability  limits  is  one  of  the  fundamental  areas 
of  combustion  science,  and  information  on  the  flammability  limits 
of  mixtures  of  fuel  and  oxidant,  with  or  without  diluents,  is  impor¬ 
tant  for  the  prevention  and  suppression  of  accidental  fires.  In  this 
study,  the  maximum  concentration  of  a  diluent  at  the  extinction 


CH4/Air/Diluent  Flames  @  (j)  =  1 .0  &  NTP 


Fig.  2.  Measured  (symbols)  and  computed  (lines)  Markstein  lengths  as  functions  of 
the  mole  fractions  of  helium,  argon,  nitrogen  and  carbon  dioxide  for  adiabatic, 
premixed  CH4/air/diluent  flames  at  ^  =  1 .0  and  NTP. 


CH4/Air/Diluent  Flames  (g  (j)  =  1 .0  &  NTP 


Fig.  3.  Measured  (symbols)  and  computed  (lines)  Markstein  numbers  as  functions 
of  the  mole  fractions  of  helium,  argon,  nitrogen  and  carbon  dioxide  for  adiabatic, 
premixed  CH4/air/diluent  flames  at  ^  =  1.0  and  NTP. 

limit  is  of  interest  because  it  is  a  measure  of  the  efficiency  of  the 
diluents  for  preventing  or  suppressing  a  fire.  The  extinguishing  dil¬ 
uent  concentration  was  determined  by  systematically  varying  the 
composition  and  ignition  energy  and  the  limiting  condition  was 
found  by  successive  experiment  trials.  A  successful  ignition  means 
that  a  flame  kernel  initially  produced  by  the  ignition  spark  can 
propagate  throughout  the  combustion  chamber.  Failure  to  ignite 
means  the  hot  gas  kernel  from  the  ignition  spark  quickly  dissipates 
and  no  flame  propagation  is  observed.  We  also  used  a  pressure 
transducer  to  record  the  pressure  history  within  the  chamber.  A 
successful  ignition  was  accompanied  by  a  rapid  increase  of  pres¬ 
sure  inside  the  chamber,  whereas  failure  to  ignite  was  accompa¬ 
nied  by  no  pressure  increase  after  the  ignition  spark. 
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Table  2 

Extinguishing  concentration  of  the  diluents  for  CH4/air  flames  at  ^  =  1 .0  and  NTP. 


Diluent 

Xd"  by 
this  work 

Xd  by  Coward  and 
Jones  [49] 

Xd  by  Liao 
et  al.  [2] 

Xd  by  Yamaoka 
and  Tsuji  [50] 

He 

38 

39 

- 

- 

Ar 

51 

48 

54 

- 

N2 

39 

37 

40 

38 

CO2 

22 

23 

29 

- 

^  Xd  is  the  diluent  concentration  in  the  reactant  mixture. 


~730  times  higher  ignition  energy  than  the  argon  mixture  if  we  as¬ 
sume  that  both  flames  have  the  same  velocity.  In  addition,  a  recent 
theoretical  study  of  Chen  and  Ju  [54]  shows  that  at  intermediate  Le¬ 
wis  numbers,  Eign  increases  exponentially  with  the  increase  of  Le. 
Therefore,  the  present  measured  limit  may  be  the  ignition  limit 
due  to  the  extremely  high  ignition  energy  required. 

Besides  the  minimum  ignition  energy,  flame/stretch  interaction 
may  also  contribute  to  the  complex  limit  behavior  of  helium.  The 
instantaneous  effect  of  stretch  on  the  flame  is  controlled  by  a  gen¬ 
eralized  stretch  factor  [55] 


Table  2  shows  the  measured  concentration  of  each  diluent  at 
extinction  for  stoichiometric  CH4/air  flames  at  NTP,  considering  he¬ 
lium,  argon,  nitrogen,  and  carbon  dioxide  as  diluents.  Xq  is  the  dil¬ 
uent  molar  concentration  in  the  reactant  mixtures.  Also  shown  in 
Table  2  are  the  flammability  limits  measured  by  Coward  and  Jones 
using  the  explosion  burette  apparatus  developed  at  US  Bureau  of 
Mines  [49],  by  Liao  et  al.  [2]  using  a  tubular  flame  burner,  and  by 
Yamaoka  and  Tsuji  [50]  using  twin  flames,  respectively.  The  pres¬ 
ent  limits  are  close  to  those  of  Coward  and  Jones,  but  systemati¬ 
cally  lower  than  those  of  Liao  et  al.  and  Yamaoka  and  Tsuji  for 
most  cases.  From  a  theoretical  point  of  view,  limits  arise  because 
mechanisms  such  as  chain-termination  reactions,  heat  losses  due 
to  conduction  or  radiation,  and  preferential-diffusion,  eventually 
dominate  the  energy-releasing  chemical  reactions  and  cause 
extinction  at  the  flammability  limit.  Flammability  limits  also  de¬ 
pend  on  the  apparatus.  The  broader  flammability  limits  using  the 
tubular  burner  method  or  the  twin  flame  method,  as  shown  in  Ta¬ 
ble  2,  may  be  because  there  is  less  heat  loss  in  the  burnt  gas  region 
than  in  the  case  of  the  spherical  flames.  In  addition,  spherically- 
expanding  flames  have  larger  stretch  rates  at  the  early  propagation 
stage,  which  may  also  explain  the  lower  limits.  The  slight  differ¬ 
ence  between  the  present  measurements  and  those  by  Coward 
and  Jones  may  be  due  to  gravity,  as  buoyant  convection  may  affect 
extinction  limits  [51,52]. 

If  the  effectiveness  of  the  diluents  are  evaluated  based  on  their 
extinguishing  concentrations,  the  ranking  is  Ar  <  He  N2  <  CO2. 
CO2  is  the  most  effective  with  the  lowest  concentration  mainly 
due  to  its  large  heat  capacity.  The  limit  behavior  of  helium  is  more 
complex.  It  is  as  effective  as  nitrogen  but  more  effective  than  ar¬ 
gon.  As  shown  in  Fig.  1,  for  a  given  diluent  concentration,  he¬ 
lium-diluted  flames  have  higher  burning  velocities  than  argon 
flames  due  to  higher  mass  and  thermal  diffusivity.  The  results  in 
Table  2,  however,  show  that  helium  is  more  effective  than  argon 
in  terms  of  the  extinguishing  concentration.  Note  the  limiting  Ar, 
N2  and  CO2  mixtures  have  a  Le  1,  while  the  limiting  helium  mix¬ 
ture  has  a  Le  4,  as  shown  in  Table  1.  When  the  helium  concentra¬ 
tion  is  above  38%,  we  observed  that  the  ignition  kernel  stopped 
propagating  very  early,  even  though  the  ignition  energy  was  signif¬ 
icantly  increased. 

To  evaluate  the  minimum  ignition  energy  required  to  ignite  a 
premixed  mixture,  the  classical  ignition  analysis  based  on  flame 
kernel  propagation  [53]  was  used  here,  which  shows  that  the  min¬ 
imum  ignition  energy  can  be  expressed  as 
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where  Cp  is  the  heat  capacity,  Tb  and  is  the  temperature  of  burned 
and  unburned  gas,  Rb  is  the  critical  radius,  a  is  the  thermal  diffusiv¬ 
ity,  and  Sl  is  the  laminar  burning  velocity.  The  quantity  Rb  is  defined 
such  that  a  flame  will  not  propagate  if  the  actual  radius  is  smaller 
than  this  critical  value.  Eign  to  be  supplied  by  the  spark  is  the  energy 
required  to  heat  the  critical  gas  volume  from  its  initial  state  to  the 
flame  temperature,  and  Eign  ~  assuming  other  parameters  are 
constants.  Therefore,  because  ane  is  approximately  nine  times  high¬ 
er  than  aAn  the  near-limit  helium  mixture  would  require  about 


where  Ze,  Le,  and  Ka  are  the  Zel’dovich  number,  the  Lewis  number, 
and  the  Karlovitz  number,  respectively.  The  flame  temperature  and 
propagating  speed  are  increased  from  the  unstretched,  adiabatic 
values  for  a  >  0,  and  decreased  for  a<0.  Since  the  outwardly-prop¬ 
agating  spherical  flames  have  positive  stretch  rate,  the  above  equa¬ 
tion  indicates  that  the  flame  temperature  and  flame  speed  are 
decreased  for  Le  >  1  mixtures  and  increased  for  Le  <  1  mixtures. 
Therefore,  for  the  near-limit  helium  mixtures  (Le  >  1 ),  stretch  sup¬ 
presses  flame  propagation.  A  competition  between  stretch  and  en¬ 
ergy  supplied  by  the  ignition  source  can  cause  the  extinction  of 
the  propagating  flame  if  the  ignition  energy  is  not  sufficient  to  drive 
the  flame  front  to  a  radius  at  which  stretch  effects  are  reduced. 

Overall,  extinction  of  the  flames  by  CO2,  N2,  and  Ar  is  likely 
caused  by  the  flame  temperature  dropping  and  radiative  heat 
losses,  and  thus  chemical  reactions  are  not  sustained.  Flame 
extinction  by  helium,  however,  is  mainly  due  to  the  high  Lewis 
number  effect.  Due  to  the  high  diffusivity  of  helium,  much  higher 
minimum  ignition  energy  is  required.  Furthermore,  stretch  sup¬ 
presses  flame  propagation  and  can  cause  flame  extinction. 


4.4.  Effects  of  radiative  emission  and  absorption 

It  is  well  known  that  radiative  heat  loss  becomes  increasingly 
important  near  flammability  limits.  Radiation,  combined  with 
stretch,  can  change  flame  dynamics  and  extinguishing  processes. 
Additionally,  recent  studies  [34]  have  shown  that  the  reabsorption 
effect  of  CO2  can  be  significant  in  some  situations  and  an  optically- 
thin  model  may  not  be  appropriate.  Motivated  by  this,  we  used  two 
radiation  models  in  the  present  COSILAB  simulations:  an  optically- 
thin  limit  model  and  an  optically-thick  model  that  solves  the  radi¬ 
ative  transfer  equation  using  DTM.  For  the  radiation  property  of 
gas  species,  gray  gas  assumption  was  used  for  the  optically-thin 
limit  model  and  a  modified  wide  band  model  was  used  for  the  opti¬ 
cally-thick  model.  In  the  following,  our  discussion  will  focus  on  the 
C02-diluted  flames  as  CO2  is  a  strong  emitter  and  absorber  com¬ 
pared  to  the  other  diluents. 

Fig.  4  compares  the  measured  and  computed  burning  velocities 
of  the  CH4/air/C02  flames  at  various  diluent  concentrations.  It  can 
be  seen  that  the  optically-thin  model  underpredicts  flame  speed. 
The  DTM  method  with  a  modified  wide  band  model  has  better 
agreement  with  the  measurements  in  the  near-limit  region.  It, 
however,  still  overpredicts  flame  speed.  The  difference  between 
the  computed  flame  speeds  using  the  two  radiation  models  in¬ 
creases  with  increasing  CO2  concentration. 

To  understand  the  effects  of  radiative  emission  and  absorption 
of  CO2  on  the  transient  flame  propagation  process,  we  chose  to  dis¬ 
cuss  a  near-limit  CH4/air  flame  with  20%  CO2  addition.  Fig.  5  shows 
the  measured  and  computed  flame  propagation  speed  relative  to 
burned  gases  as  a  function  of  radius.  The  computed  burned  gas 
flame  speeds  were  obtained  with  an  optically-thin  model  (dashed 
line)  and  with  the  DTM  model  (solid  line).  It  can  be  seen  that  radi¬ 
ation  does  not  play  an  important  role  during  the  early  stage  of 
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Fig.  4.  Measured  (symbols)  and  computed  (lines)  laminar  burning  velocities  as  a 
function  of  the  mole  fraction  of  carbon  dioxide  for  premixed  CH4/air/C02  flames  at 
^  =  1.0  and  NTP.  Solid  line:  DTM  method  and  a  modified  wide  band  model;  dashed 
line:  optically-thin  limit  model. 


Fig.  5.  The  measured  (symbols)  and  computed  (lines)  flame  propagation  speed 
relative  to  burned  gases  as  a  function  of  flame  radius  for  the  CH4/air/20%  CO2  flame 
at  ^  =  1.0  and  NTP.  Solid  line:  DTM  method  and  a  modified  wide  band  model; 
dashed  line:  optically-thin  limit  model. 


flame  propagation  when  r  is  small,  e.g.,  r<1.5  cm.  Radiation 
becomes  increasingly  important  when  r  becomes  large.  After  the 
initial  flame  propagation  stage,  the  flame  speed  predicted  by  the 
optically-thin  limit  model  decreases  due  to  heat  losses  to  the  cold 
ambient  gases.  On  the  contrary,  the  flame  speed  predicted  by  the 
optically-thick  model  increases  with  radius.  This  is  mainly  because 
the  unburned  gases  absorb  heat  from  the  hot  zone,  which  results  in 
elevated  temperature  of  the  unburned  gases  and  thus  increases 
flame  propagation  speed.  The  optically-thick  model  still  overpre¬ 
dicts  flame  speed,  indicating  that  a  more  sophisticated  radiation 
property  model  such  as  a  narrow  band  model  may  be  needed. 

The  observed  variation  of  flame  speed  due  to  radiative  emission 
or  absorption  also  raises  a  concern  on  the  linear  extrapolation 
method  which  has  been  widely  used  to  extract  unstretched  burn¬ 
ing  velocities  from  spherically-expanding  flames.  Strictly  speaking, 
the  linear  relation  (Eq.  (4))  between  stretched  flame  speed  and 
stretch  rate  is  based  on  the  assumption  of  small  stretch  and  adia¬ 
batic  flames.  For  radiating  near-limit  flames,  however,  flame  speed 
variation  due  to  radiation  does  not  have  a  linear  dependence  on 
flame  stretch.  Therefore,  for  such  flames,  a  non-linear  extrapola¬ 
tion  method  that  includes  the  effect  of  radiation  may  be  more 
appropriate. 


Radius  (cm) 


Fig.  6.  The  computed  temperature  (top)  and  radiative  power  (bottom)  as  a  function 
of  radius  at  various  times  for  the  spherically  expanding  CH4/air/20%  CO2  flame  at  (j) 
=  1.0  and  NTP  using  an  optically-thin  radiation  model.  Ignition  occurs  at  t  =  0. 
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Fig.  7.  The  computed  temperature  (top)  and  radiative  power  (bottom)  as  a  function 
of  radius  at  various  times  for  the  spherically  expanding  CH4/air/20%  CO2  flame  at  (j) 
=  1.0  and  NTP  using  an  optically-thick  model  based  on  the  DTM  method  and  a 
modified  wide  band  model.  Ignition  occurs  at  t  =  0. 

Figs.  6  and  7  show  the  computed  profiles  of  temperature  and 
radiative  power  after  ignition  for  the  same  near-limit  flame  as 
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discussed  above  (CH4/air/20%  CO2)  using  the  optically-thin  radia¬ 
tion  model  (Fig.  6)  and  the  DTM  model  (Fig.  7)  respectively.  As 
shown  in  Fig.  6,  beyond  the  initial  ignition  stage,  the  temperature 
of  burned  gases  continues  to  decrease  as  the  flame  propagates  out¬ 
wardly  due  to  radiative  heat  loss  to  the  cold  ambient  gases.  The 
radiative  power  exhibits  similar  trend  as  the  temperature.  The 
peak  temperature  of  the  flame  zone  is  around  1650  K,  with  only 
slight  decrease  (3-5  K)  for  each  of  the  next  time  steps.  Similar 
trends  were  observed  in  Fig.  7  using  the  DTM  method  and  a  mod- 
ifled  wide  band  model:  the  temperature  and  radiative  power  of 
burned  gases  decrease  as  flame  propagates  outwardly.  However, 
beyond  the  ignition  stage  the  peak  temperature  of  the  reaction 
zone  continues  to  increase  as  flame  propagates  outwardly.  The 
radiative  power  profiles  show  that  a  significant  amount  of  radia¬ 
tion  heat  loss  from  the  hot  products  is  reabsorbed  by  the  cold 
ambient  gases,  resulting  temperature  increase  of  the  unburned 
gases.  The  peak  radiative  power  at  the  flame  zone  remains  around 
3.6  X  10®W/m^  during  flame  propagation,  with  only  slight  de¬ 
crease  for  each  of  the  next  time  steps.  The  radiative  power  ab¬ 
sorbed  by  the  cold  unburned  gases,  however,  continues  to 
increase  with  increasing  r,  due  to  the  increasing  volume  of  the 
hot  burned  gases.  When  the  flame  propagates  to  r  =  4.5  cm,  the 
peak  radiation  power  absorbed  by  the  cold  gases  is  as  large  as 
the  peak  power  emitted  from  the  hot  gases. 

4.5.  Chemical  effect  of  CO2  in  flame  suppression 

The  inert  gases  considered  here  belong  to  chemically  passive  or 
thermal-quenching  types  of  flame  suppressants,  for  which  flame 
temperature,  flame  speed  and  reaction  rate  are  reduced  mainly 
due  to  dilution  and  thermal  effects.  This  inhibition  mechanism  is 
different  from  chemically  active  types  of  suppressants  such  as 
Halons  and  Halon  replacements.  Active  flame  suppressants  reduce 
reaction  rate  mainly  by  reducing  the  concentration  of  chain-carry¬ 
ing  radicals  (such  as  H  and  OH)  in  the  reaction  zone.  For  example, 
Halon  1301  (CF3Br)  consumes  H-atoms  in  the  reaction  zone 
through  the  reactions  of  CFsBr  +  H  =  CF3  +  HBr,  H  +  HBr  =  H2  +  Br 
[56]. 

While  N2,  Ar  and  He  suppress  a  flame  mainly  by  thermal  effect, 
flame  inhibition  by  CO2  involves  both  physical  and  chemical  ef¬ 
fects,  which  have  been  studied  by  several  researchers  [57-60]. 
These  studies  show  that  CO2  can  be  converted  to  CO  through  the 
reaction  CO  +  OH  =  CO2  +  H,  which  is  the  dominant  reaction  path¬ 
way  for  the  chemical  participation  of  CO2  in  hydrocarbon  flames. 
The  competition  of  CO2  for  H  radicals  through  this  reaction  with 
the  most  important  chain-branching  reaction  H  +  O2  =  0  +  OH  can 
play  a  chemical  inhibiting  role  that  reduces  the  concentrations  of 
important  radicals,  such  as  0,  H,  and  OH,  which  subsequently  re¬ 
duces  reaction  rate  and  burning  speed.  Liu  et  al.  [59]  investigated 
the  chemical  effect  of  CO2  on  premixed  CH4  and  H2  flames  numer¬ 
ically.  Their  results  show  that  the  chemical  effect  of  CO2  signifi¬ 
cantly  reduces  the  burning  velocity  and  the  relative  importance 
of  the  chemical  effect  of  CO2  increases  with  increasing  CO2 
concentration. 

To  understand  the  importance  of  chemical  effect  of  CO2  in  the 
overall  flame  suppression  process,  we  did  a  numerical  analysis  to 
quantify  the  pure  chemical  effect  of  CO2  in  flame  extinction.  The 
approach  is  similar  to  the  work  of  Liu  et  al.  [59]  and  Park  et  al. 
[58,60].  This  was  done  using  the  1-D  planar  freely-propagating 
flame  code  PREMIX  [16].  The  thermodynamics  and  transport  data 
were  taken  from  CHEMKIN  library  [41,42].  The  PRMEIX  calcula¬ 
tions  used  the  same  reaction  mechanism,  GRl-Mech  3.0  [43],  as 
the  transient  COSILAB  simulations.  The  strategy  was  to  introduce 
an  artificial  species,  named  X,  which  has  the  same  thermochemical, 
transport,  and  radiation  properties  as  CO2.  However,  this  species 
does  not  participate  in  chemical  reactions  and  is  chemically  inert. 


Fig.  8.  The  computed  laminar  burning  velocity  and  burned  gas  temperature  as  a 
function  of  CO2  mole  faction  for  the  adiabatic  CH4/air/C02  flames  at  ^  =  1 .0  and  NTP 
using  PREMIX.  Solid  line:  results  for  real  CO2:  dashed  line:  results  for  artiflcial  CO2 
which  does  not  participate  in  chemical  reactions. 


Fig.  9.  The  computed  profiles  of  temperature  and  H,  0,  and  OH  radical  concentra¬ 
tions  for  the  adiabatic  CH4/air/20%  CO2  flame  at  0  =  1.0  and  NTP  using  PREMIX.  Solid 
lines:  results  for  real  CO2:  dashed  lines:  results  for  artiflcial  CO2  which  does  not 
participate  in  chemical  reactions. 

Two  solutions  were  obtained  for  each  C02-diluted  flame:  one  for 
real  CO2  which  participates  in  chemical  reactions,  the  other  for 
the  artiflcial  species  X  excluding  chemical  effect.  In  addition,  the 
species  X  has  the  same  third-body  collision  efficiency  as  CO2  in 
all  relevant  reactions. 

Eig.  8  shows  the  computed  laminar  burning  velocity  and  the 
temperature  of  burned  gases  of  C02-diluted  CH4/air  flames  at 
NTP  and  (/>  =  1.0  as  a  function  of  the  mole  fraction  of  CO2  using 
the  method  described  above.  The  solid  line  represents  the  results 
for  real  CO2  and  the  dash  line  represents  the  results  for  artiflcial 
CO2.  In  these  calculations,  radiation  heat  loss  and  reabsorption 
were  not  considered.  As  can  be  seen  in  Eig.  8,  the  chemical  effect 
causes  significant  reduction  in  burning  velocity  and  the  degree  of 
flame  speed  reduction  increases  with  increasing  CO2  concentra¬ 
tion.  Near  the  limit  condition,  e.g.,  20%  CO2  dilution  by  volume, 
the  burning  velocities  are  7.1  cm/s  and  9.9  cm/s  for  real  CO2  and 
artiflcial  CO2  respectively.  This  means  chemical  effect  causes  about 
30%  reduction  in  flame  speed.  The  difference  in  flame  temperature, 
however,  is  not  significant,  about  30-40  K  for  all  CO2  mole 
fractions. 
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Fig.  10.  Net  reaction  rates  of  stoichiometric  CH4/air/20%  CO2  flame  at  NTP. 


We  then  plotted  the  computed  temperature  profile  and  the  pro¬ 
files  of  H,  0,  and  OH  radical  concentrations  for  the  near-limit  CH4/ 
air/20%  CO2  flame  (Fig.  9).  The  chemical  effect  of  CO2  reduces  H, 
OH,  and  0  radical  concentrations.  In  particular,  the  maximum  OH 
radical  concentration  for  the  artificial  CO2  case  is  almost  twice  as 
for  the  real  CO2  case.  The  radical  reductions  are  mainly  due  to 
the  reaction  CO  +  OH  =  CO2  +  H,  for  which  the  reaction  rate  is  the 
second  largest  among  all  elementary  reactions  and  the  first  largest 
is  the  chain-branching  reaction  H  +  O2  =  0  +  OH,  as  shown  in 
Fig.  10.  The  competition  of  CO2  for  H  radical  through  CO  +  OH  = 
CO2  +  H  with  H  +  O2  =  0  +  OH  reduces  H,  OH,  and  0  radical  concen¬ 
trations,  leading  to  lower  reaction  rate  and  flame  speed. 

5.  Conclusions 

An  experimental  and  numerical  study  was  carried  out  to  inves¬ 
tigate  the  extinction  behavior  of  premixed  methane-air  flames  by 
various  diluents  including  helium,  argon,  nitrogen  and  carbon 
dioxide  in  microgravity.  The  experiments  were  conducted  using  a 
short-drop  free-fall  facility  that  images  flame  propagation  in  a 
dropping  chamber  with  stationary  shadowgraph  system.  This  en¬ 
ables  accurate  determination  of  the  near-limit  burning  velocity. 
The  maximum  concentration  of  a  specific  diluent  at  the  extinction 
limit  was  determined  by  systematically  varying  the  composition 
and  ignition  energy  and  finding  the  limiting  condition  through  suc¬ 
cessive  experiment  trials.  The  corresponding  freely-propagating, 
planar  flames  were  simulated  using  PREMIX.  The  spherically- 
expanding  flames  were  simulated  using  the  Spherical  Flames  Mod¬ 
ule  of  COSILAB  considering  detailed  radiation  models  in  spherical 
geometry.  The  following  conclusions  can  be  drawn  based  on  the 
present  experiments  and  computations: 

1.  Helium  exhibits  more  complex  limit  behavior  than  the  other 
diluents.  Extinction  of  the  flames  by  CO2,  N2,  and  Ar  is  likely 
caused  by  the  flame  temperature  dropping  and  radiative  heat 
losses,  and  thus  chemical  reactions  are  not  sustained.  Elame 
extinction  by  helium,  however,  is  mainly  due  to  the  high  Lewis 
number  effect.  Due  to  the  high  thermal  diffusivity,  much  higher 
ignition  energy  is  required  for  helium  mixtures  than  for  the  oth¬ 
ers.  In  addition,  for  spherically  expanding  helium-diluted 
flames,  stretch  suppresses  flame  propagation  and  can  cause 
flame  extinction. 

2.  Both  an  optically-thin  model  and  an  optically-thick  model  were 
employed  in  the  spherically-expanding  flame  simulations  using 
COSILAB  to  understand  the  effects  of  radiative  emission  and 
absorption  on  flame  propagation  and  extinction.  The  results 
show  that  the  optically-thin  model  underpredicts  flame  speed 
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and  the  optically-thick  model  based  on  the  DTM  method  and 
a  modified  wide  band  model  has  better  agreement  with  mea¬ 
surements  in  the  near-limit  region.  The  simulation  results  show 
that  a  significant  amount  of  radiation  heat  loss  from  the  hot 
products  is  reabsorbed  by  the  cold  ambient  gases,  resulting  in 
temperature  increase  in  the  unburned  gas  and  the  pre-heat 
zone.  The  peak  radiative  power  due  to  absorption  increases 
with  time.  The  optically-thick  model,  however,  still  overpre¬ 
dicts  flame  speed,  indicating  that  a  more  sophisticated  radiation 
property  model  such  as  a  narrow  band  model  may  be  needed.  In 
addition,  a  non-linear  extrapolation  method  that  includes  the 
effect  of  radiation  may  be  more  appropriate  for  C02-diluted 
near-limit  flames. 

3.  While  N2,  Ar  and  He  suppress  a  flame  mainly  by  thermal  effect, 
flame  inhibition  by  CO2  involves  both  physical  and  chemical 
effects.  A  numerical  analysis  shows  that  the  pure  chemical 
effect  of  CO2  causes  about  30%  reduction  in  flame  speed  for 
the  near-limit  CH4/air/20%  CO2  flame.  The  reduction  is  because 
CO2  participates  in  the  reaction  CO2  +  H  =  CO  +  OH,  which  com¬ 
petes  for  H  radicals  with  the  most  important  chain-branching 
reaction  H  +  O2  =  0  +  OH  and  thus  reduces  the  concentrations 
of  important  radicals,  such  as  0,  H,  and  OH.  This  effect  becomes 
more  important  when  the  CO2  concentration  becomes  large. 
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